Kondo Physics in the Single Electron Transistor with ac Driving 
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Using a time-dependent Anderson Hamiltonian, a quantum dot with an ac voltage applied to a nearby gate is investigated. 
A rich dependence of the linear response conductance on the external frequency and driving amplitude is demonstrated. At low 
frequencies the ac potential produces sidebands of the Kondo peak in the spectral density of the dot, resulting in a logarithmic 
decrease in conductance over several decades of frequency. At intermediate frequencies, the conductance of the dot displays 
an oscillatory behavior due to the appearance of Kondo resonances of the satellites of the dot level. At high frequencies, the 
conductance of the dot can vary rapidly due to the interplay between photon-assisted tunneling and the Kondo resonance. 

PACS numbers: 72.15.Qm, 85.30.Vw, 73.50.Mx 
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It has been predicted that, at low temperatures, trans- 
port through a quantum dot should be governed by the 
same many-body phenomenon that enhances the resistiv- 
ity of a metal containing magnetic impurities - namely 
the Kondo effect Bj . The recent observation of the Kondo 
effect by Goldhaber et al. 0] in a quantum dot operat- 
ing as a single-electron transistor (SET) has fully verified 
these predictions. In contrast to bulk metals, where the 
Kondo effect corresponds to the screening of the free spins 
of a large number of magnetic impurities, there is only 
one free spin in the quantum-dot experiment. Moreover, 
a combination of bias and gate voltages allow the Kondo 
regime, mixed- valence regime, and empty-site regime all 
to be studied for the same quantum dot, both in and out 
of equilibrium 

Here we consider another opportunity presented by the 
observation of the Kondo effect in a quantum dot that 
is not available in bulk metals - the application of an 
unscreened ac potential. There is already a large liter- 
ature concerning the experimental application of time- 
dependent fields to quantum dots ||. For a dot acting 
as a Kondo system, the ac voltage can be used to peri- 
odically modify the Kondo temperature or to alternate 
between the Kondo and mixed-valence regimes. Thus it 
is natural to ask what additional phenomena occur in a 
driven system which in steady state is dominated by the 
Kondo effect ||Jj|. Our results indicate a rich range of 
behavior with increasing ac frequency, from sidebands of 
the Kondo peak at low ac frequencies, to conductance 
oscillations at intermediate frequencies, and finally to 



photon-assisted tunneling at high ac frequencies. 

The system of interest is a semiconductor quantum 
dot, as pictured schematically in Fig. [|. An electron can 
be constrained between two reservoirs by tunneling bar- 
riers leading to a virtual electronic level within the dot at 
energy ~£dot (measured from the Fermi level) and width 
~2Pdot @- We assume that both the charging energy 
e 2 /C and the level spacing in the dot are much larger 
than Tdot, so the dot will operate as a SET ||. In this 
work, we consider only the linear-response conductance 
between the two reservoirs. However, we will allow an 
oscillating gate voltage V s (t) — Vq + V^ccosfii of arbi- 
trary (angular) frequency Q and arbitrary amplitude V ac , 
which modulates the virtual- level energy £dot(0- 

Such a system may be described by a constrained (U = 
oo) Anderson Hamiltonian 



/dot(i) n a + 

kcr 



^kcr^ko 



H.c.) . (1) 



Here c\ creates an electron of spin a in the quantum dot, 
while n a is the corresponding number operator; c\ a cre- 
ates a corresponding reservoir electron; k is shorthand 
for all other quantum numbers of the reservoir electrons, 
including the designation of left or right reservoir, while 
14 is the tunneling matrix element through the appropri- 
ate barrier. Because the charging energy to add a second 
electron, U = e 2 /C, is assumed large, the Fock space in 
which the Hamiltonian ([!]) operates is restricted to those 
elements with zero or one electron in the dot. 
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FIG. 1. Schematic picture of the quantum dot SET. 

At low temperatures, the Anderson Hamiltonian (|l|) 
gives rise to the Kondo effect when the level energy edot 
lies below the Fermi energy. In this regime, a single elec- 
tron occupies the dot which, in effect, turns the dot into 
a magnetic impurity with a free spin. The temperature 
required to observe the Kondo effect in linear response 
is of order Tk ~ £>exp(— 7r|edot|/rdot), where D is the 
energy difference between the Fermi level and the bot- 
tom of the band of states. For the temperature range 
that is likely to be experimentally accessible in a SET, 
T ~ Tk or higher, there exists a well tested and reliable 
approximation known as the non-crossing approximation 
(NCA) 0. The NCA has been formally generalized to 
the full time-dependent nonequilibrium case and an 
exact method for the (numerical) solution implemented 
H . The time-dependent NCA has been applied to Kondo 
physics in charge transfer in hyperthermal ion scatter- 
ing from metallic surfaces and to energy transfer 
and stimulated desorption at metallic surfaces 11 1. An 
independent formulation Q has been applied to quan- 
tum dots, although the high frequency expansion used 
there appears limited. Here we present the exact time- 
dependent NCA solution for a quantum dot over the full 
range of applied frequencies. 

The time-dependent electronic structure of the dot can 
be characterized by the time-dependent spectral density 



Pdot{e,t) 



(*-|t)}) (2) 



evaluated in the restricted Fock space. For the equi- 
librium Kondo system, /Odot(e) is time independent, and 
looks like the graph in the schematic in Fig. [|. Roughly 
speaking, pdot(e) consists of a broad peak of width 
^2F(jot at the level position edot and a sharp Kondo peak 
of width ~Xk near the Fermi level. We will refer to these 
features as the virtual-level peak and the Kondo peak, re- 
spectively. In the steady-state case, the linear-response 
conductance G through a dot symmetrically coupled to 
two reservoirs is given by fll2|] 



G = 



dot 



h 2 



de pdot(e) 



dm 

de 



(3) 



where /(e) is the Fermi function. The formula (^|) will 
still be valid in the case where the gate voltage is time 
dependent if G is the time-averaged conductance and 



Pdot(e) is replaced by the time-averaged spectral density 
(Pdot(e,i)) F° r a given system, this average will 

depend on the driving amplitude V ac and frequency O. 
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FIG. 2. The spectral density {p^ ot (e,t)) vs. energy e for a 
quantum dot with level energy edot(i) = —5 + 4 cos fit and 
T — 0.005. The non-driven case is also shown in the final 
panel. Throughout this letter, energies are in units of Tdot- 

In Fig. H we show the calculated (pdot(e, t)) as a func- 
tion of energy e for a level with energy edot(i) = fdot + 
e ac cosfii at several different frequencies fl. The corre- 
sponding conductance is shown by the curve labeled dot 
A, T = 0.005 in Fig. § For the lowest fi, the response of 
the system is relatively adiabatic and the displayed spec- 
tral function resembles the spectral function that would 
have resulted if the system had been in perfect equilib- 
rium for all the dot level positions over a period of oscil- 
lation of edot {t) ■ The two broad peaks are the influence 
of the virtual level peaks at the two stationary points of 
this oscillation (here at e = — 1 and e = — 9). As the fre- 
quency f2 is increased, marked nonadiabatic effects result, 
the most obvious being the appearance of multiple satel- 
lites around the Kondo resonance These sidebands 
appear at energies equal to Ti times multiples of the driv- 
ing frequency fl |fL5| . As the frequency Cl is increased, 
spectral weight is transferred from the main Kondo peak 
to these satellites. As the conductance is dominated by 
(Pdot(eji)) at the Fermi energy, this causes the slow log- 
arithmic falloff of the conductance over two decades of 
frequency, as shown in Fig. ||[ 

As Ml becomes larger than Tdot, inspection of Fig. || 
shows that broad satellites also appear at energy separa- 
tions nhfl around the average virtual- level position edot- 
These satellites of the virtual level are the analogues of 
those predicted in the noninteracting case p~3| , which de- 
crease in magnitude as the order n of the Bessel function 
J n . Here, however, the virtual-level satellites have their 
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own Kondo peaks; each of the latter gets strong when 
the corresponding virtual-level satellite reaches a posi- 
tion a little below the Fermi level, and then disappears 
as the broad satellite crosses the Fermi level. This effect 
produces the oscillations in the conductance that are ev- 
ident in the lower curves in Fig. ||. These oscillations are 
very different from those that would occur in a noninter- 
acting (U = 0) case: due to the Kondo peaks they are 
substantially stronger, their maxima occur at different 
frequencies and their magnitudes are temperature depen- 
dent. As the last virtual- level satellite crosses the Fermi 
level, Ml = |edot|, the dot level energy begins to vary 
too fast for the system to respond and the average spec- 
tral function approaches (exactly as f2 — > oo) the equi- 
librium spectral function for a dot level centered at the 
average position edot- For the parameters of Fig. |[ the 
high frequency region is uninteresting, because the tem- 
perature is far above the Kondo temperature. Therefore, 
the conductance shows little temperature or frequency 
dependence at these high frequencies. 




FIG. 3. Conductance of two different quantum dots, each 
at two different temperatures: Dot A, edot(t) = — 5 + 4cosf2t; 
Dot B, £dot(i) = —2.5 + 2 cos fit. The curves at the 
high f2 end for dot B (marked "PAT") are from our pho- 
ton-assisted-tunneling model, while the exact high frequency 
asymptotes for dot B are shown as short horizontal lines ex- 
tending from the right vertical axis. The inset shows the 
spectral density (pdot(e, t)) of dot B around the Fermi level, 
at T = 0.005, for large frequencies, from f2 = 4.8 (lowest 
curve), through 5.5, 6.1, 6.8 to f2 = 14 (topmost curve). 

The situation is quite different for the system (dot B) 
displayed in the upper two curves in Fig. |^, which dis- 
plays a strong Kondo effect when the dot level is held at 
its average energy edot- In this case the fl — > oo conduc- 
tance is strongly enhanced by the Kondo effect, and is 
consequently temperature dependent as well. Note that 
the conductance falls off significantly from its asymptotic, 



Q, — > oo, value for frequencies still much larger than ei- 
ther the depth of the level edot or its width. This effect 
is due to a rapid decline of the amplitude of the Kondo 
peak in the spectral density, as illustrated in the inset 
of Fig. §. We propose the following explanation for this 
phenomenon. The energy Ml excites the dot, produc- 
ing satellites M,ff3[ of the virtual level peak at energies 
edot i nhil, which, for Ml 3> Tdot have strength roughly 
given by [J ra (e ac /?ifi)] 2 as in the U = case (see Fig. || 
and the previous discussion). For large hfl, only the two 
n = 1 satellites have any significant strength, and the 
higher lies above the Fermi level, allowing an electron on 
the dot to decay at the rate (l/?i)rdot(edot + hCl). The 
overall electron decay probability per unit time Tdccay/^ 
due to this photon-assisted-tunneling mechanism (PAT) 
is therefore given by 

Tdccay ~ [</i(<W?i^)] 2 rdot(edot + MI). (4) 

The above rate carries with it an energy uncertainty, 
which we speculate has roughly the same effect on the 
Kondo peak as the energy smearing due to a finite tem- 
perature. We can test this conjecture by calculating the 
equilibrium conductance at an effective temperature T e g 
given by T c g = T + Tdocay The results of such a cal- 
culation are shown in Fig. |^ (PAT curves), where they 
compare very favorably with our results for the conduc- 
tance in the ac-driven system. 

Returning to the behavior at low frequencies, we find 
that it can be best understood in terms of the Kondo 
Hamiltonian, which, with respect to properties near the 
Fermi level, is equivalent to the Anderson Hamiltonian 
(|l|) in the extreme Kondo region — edot S> Tdot 01- In 
this limit the dot can be replaced simply by a dynami- 
cal Heisenberg spin S (S 2 = |), which scatters electrons 
both within and between reservoirs. The Kondo Hamil- 
tonian corresponding to the Anderson model |l]) is 

c L c fc'<7'> ( 5 ) 

kk' oo' 

where the components of a are the Pauli spin matri- 
ces. For near Fermi level properties we can suppress 
the detailed k dependence of J and V and introduce a 
large energy cutoff D m, 

in which case the relation- 
ship between the Kondo and Anderson Hamiltonians Jl]]] 
is J(i) = |V^ 2 /edot(£)| for our U = oo case. If we let 
wieads(e) / "h be the total rate at which lead electrons of 
energy e undergo intralead and interlead scattering by 
the dot, then wi ea ds(e) will have a Kondo peak for e 
near the Fermi level. Furthermore, if J is modulated 
as J(t) = (J) (1 + a cos ilt), then an electron scattered by 
the dot will be able to absorb or emit multiple quanta 
of energy hil, leading to satellites of the Kondo peak in 
(pdot(e,i)) through the exact Anderson model relation 
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FIG. 4. Spectral density (pdot(e, t)} times 10 3 in the Kondo 
model and NCA for T = 0.02 and (J)p = 0.023 (e do t = -7). 
For the non-driven case (left panel) we also show the compa- 
rable result from summing all the leading logarithmic terms 
(Abrikosov, Ref. jl9|), as well as that obtained to order J 2 . 
The energy dependence of (J) (Ref. JlTj]) has been included 
to order J 2 in all the Kondo Hamiltonian curves. 



wieads(e)- This then reflects back on 

Wlcads(e) = r d ot(e)(pdot(e,t)} /pleads (e), 



(6) 



where pieads(e) is the state density per spin in the leads. 

The above can be illustrated explicitly using pertur- 
bation theory in J. Keeping all terms of order J 2 and 
logarithmic terms to order J 3 , we find, using a nonequi- 
librium version of Abrikosov's pseudofermion technique 
111, that 



Mdeads(e) = 2ir(J 2 )p 



l + 3(J)p a n g(e + nhQ) 



n=-l 



(7) 



where p = pi ca ds(0), a = 1, a ±1 = a 2 / (2 + a 2 ), (J 2 ) 
i 

" V ,1-2/(0 



(l + ±a 2 )(J) 2 , and 



de' 



In 



— D 



(8) 



the last limit being approached when T -C |e|. In Fig. |] 
we compare this prediction with the full NCA theory. 
Although we are not strictly in the parameter region 
where the J 3 theory is quantitatively valid, the quali- 
tative agreement is quite satisfactory. 

The present results indicate rich behavior when an 
external ac potential is applied to a quantum dot in 
the regime where the conductance is dominated by the 
Kondo effect. While the time-dependent NCA method 
employed spans the full range of applied frequency, some 
additional insight has been gained into the behavior both 
at very low and very high frequencies. At low frequencies 
a time-dependent Kondo model helps explain the ampli- 
tudes of sidebands of the Kondo peak in the spectral 
density of the dot. At high frequencies, a cutoff of the 



Kondo peak due to photon-assisted tunneling processes 
accounts for the reduction of conductance. We hope that 
our work will inspire experimental investigation of these 
phenomena and other ramifications of ac driving applied 
to Kondo systems. 
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